Disentanglement Induced Magnetoresistance in Organic Semiconductors 
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By considering explicitly the contribution of carrier's spin in various incoherent charge transfer 
(CT) process, we propose a fundamental mechanism for magnetoresistance in organic semiconduc- 
tors. It is proven that the disentanglement between the carrier's spin and the local environment, 
accompanying the incoherent CT, contributes a reduction factor to the intermolecular CT rate. The 
entanglement is altered by the applied magnetic field, which leads to the magnetoresistance. Based 
on a hyperfine interaction model, we show the CT rate exhibits a general feature of a Lorenzian- 
shape saturation at large fields, together with an ultrasmall-field component. Furthermore, the 
ultrasmall-field component is shown to be a result of the coherent motion of carrier's spin, while the 
saturation behavior survives the loss of coherence. The magnetic-field dependence, as well as the 
obtained isotope effect, are in satisfying agreement with the experimental results [T.D. Nguyen et 
al., Nature Materials 9, 345 (2010)]. The mechanism also explains the diversity of magnetoresistance 
observed in various organic materials. 

PACS numbers: 72.20.My, 72.80.Le, 31.30.Gs 
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Recently, a series of experiments have revealed the 
surprising long-range quantum coherence at picosecond 
timescales in the excitation energy transfer among the 
light-absorbing molecules in photosynthetic proteins 
It is further proposed that the the light-harvesting ef- 
ficiency can be enhanced by the accompanying quan- 
tum entanglement 0, whose essential role is also sup- 
ported by a theoretical analysis concluding that the en- 
tanglement could take effect for nonequilibrium systems 
even at room temperature Q ■ These discoveries indicate 
that, the processes related to quantum coherence, such as 
quantum entanglement and disentanglement, which are 
previously thought to be fragile, could lead to new phys- 
ical consequences in "natural" and molecular systems. 

As a promising branch of molecular electronics, the 
organic thin-film devices have raised much attention due 
to their low cost, flexibility and ease of fabrication [4|. 
The microscopic process of charge transport and recom- 
bination is thought to be a series of charge transfer (CT) 
between molecules, or more abstractly, transport sites, 
under different conditions. For example, a single charge 
carrier on a specific site can hop onto a neighboring va- 
cant site. Moreover, if the destination site is already 
occupied by an oppositely charged carrier, a singlet or 
triplet exciton will form. More complicated processes 
are expected such as the CT to and off a site with a 
long-lived triplet exciton, namely the triplet-polaron in- 
teraction. Due to both the amorphous morphology and 
the self-trapping effect brought about by the electron- 
lattice interaction, most of the CT processes are thought 
to be incoherent, and are well described by the CT rate. 
Various theories for the CT rate have been proposed 
for different situations, such as the Miller- Abrahams for- 
mula Q and the Marcus formula These theories, in 
combination with the Gaussian disorder model, capture 
the electric-field and temperature dependence of mobil- 
ity in a certain range successfully Q. The discovery of 



the magnetoresistance (MR), together with the magneto- 
electroluminescence (MEL) [9|,ll0( in the organic thin film 
devices raises a new topic on the role of the carrier's spin 
in the CT processes. Novel physics is expected, which 
differs from that of the inorganic counterparts due to 
the coherence of the carrier's spin and the incoherence of 
charge transport in such systems 11, 12] ■ Based on the 



experiments under different operating conditions, a num- 
ber of quantities are claimed to be modulated by mag- 
netic field [HI], including the charge carrier mobility flij j. 
the inter-system crossing rate between singlet and triplet 
electron- hole pairs [l5j |. and the recombination mobility 
[lij ]. It should be noted that, despite the diversity of 
the underlying processes, both the MR and MEL show 
surprising universality and similarity in a wide range of 
materials [ItJ and temperatures [ijj], which implies that 
a general understanding from the basic characteristics of 
the underlying CT processes is required. 

In this Letter, we study explicitly the motion of the 
carrier's spin in the incoherent CT processes of or- 
ganic semiconductors. We point out that a disentangle- 
ment process will accompany the incoherent CT, through 
which the CT rate is related to the entanglement between 
the carrier's spin and the local environment. This entan- 
glement can be modulated by the applied magnetic field, 
which leads to the MR, which is a fundamental mecha- 
nism beyond the common "spin-selection" rules. By con- 
sidering a hyperfine-interaction model and the hopping 
of a single carrier from a specific site to its neighboring 
vacant site, we calculated the magnetic-field dependence 
of the hopping rate, which shows a general feature of a 
Lorenzian-shape saturation at large fields, together with 
an ultrasmall-field component, which is observed both for 
MR and MEL [l7l[l9| . Furthermore, the ultrasmall-field 
component is shown to be a result of the spin's coherent 
motion while the saturation behavior survives the loss 
of spin coherence. The isotope effect obtained from this 
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mechanism by replacing the spin-i protons with the spin- 
1 deuterons is in satisfying agreement with experimental 
results. 

Firstly, we illustrate how the motion of the carrier's 
spin could influence the incoherent CT rate. A two-site 
model with site i and j is taken for convenience, and the 
carrier is initially localized on site i. Once on the site, 
the carrier's spin will inevitably interact with some local 
degree of freedom, which consists of the local environ- 
ment of spin (LES). This mutual interaction will bring 
the composite system of both the carrier's spin and LES 
into entangled states. The entanglement is present both 
for the case of dynamic evolution, and for the case of 
the composite system relaxing to thermal equilibrium. 
Concurrent with the interaction concerning the carrier's 
spin is the interaction between the carrier and the lattice 
vibrations, which offers the energy for the CT to site j 
with different electronic energy. The whole CT process is 
incoherent due to the small relaxation time of the lattice 
vibrations which are on the order of picoseconds [20[ . In 
terms of the density operator p of the two-site system, 
the incoherentness is reflected by the quick vanishing of 
the off-diagonal matrix elements pij (= (i \ p \j)). The in- 
coherent CT process also have substantial influence on 
the final state of the carrier's spin after the CT. As the 
interactions between the carrier and the lattice vibrations 
are not spin-specific, the coherence within the subspace 
of the the carrier's spin should remain intact to these in- 
teractions leading to fast decoherence. However, as the 
coherence between the site i and site j is lost after the 
incoherent CT, the entanglement between the carrier's 
spin (now on site j) and the LES on site i, as a sign of 
the nonlocal quantum coherence, is expected to be termi- 
nated [2~fl | . The above two requirements can be achieved 
by an adiabatic elimination procedure. If the initial state 
of the composite system for the CT is designated as pi, 
the density matrix for the carrier's spin/the LES after the 
CT is p s / e = tr e / s {pi}, where tr e / s denotes the partial 
trace over the degrees of freedom of the LES/the carrier's 
spin. The final unentangled state pf can be written as 
a direct product of the states p s and p e . By a gener- 
alization of the Fermi golden rule, the contribution to 
the CT rate from this disentanglement process can be 
incorporated in a quantity 77, which reads 



77 = tr {piPf} = tr {pi{p s <£> p e )} • 



(1) 



To show the physical meaning of 77 more clearly, we sup- 
pose the initial state pt is a pure state, which can be 
written as pi = \4>i) (</>i|, with 
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where |e+) and |e_) arc the states of the LES of site i, 
that are correlated with the states of the carrier's spin 
with magnetic quantum number m 



c + and c_ are the corresponding superposition coeffi- 
cients. In this case, 77 can be simplified to 



77 = 1 - 3|c_ 
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which is directly related to the degree of entanglement 
between the carrier's spin and the LES on site i in the 
initial state . If the contributions from other parts are 
denoted by i/q, the complete form of CT rate can be writ- 
ten as ^o?7. It should be noted that, besides the details 
of the electronic states and the interaction with lattice 
vibrations, the value of vq also depends on whether site 
j has already been occupied by a carrier, as well as the 
spin states, which has been considered by several papers 
[T^H^I- Here we focus on the influence of 77, originating 
from the disentanglement between the carrier's spin and 
the LES of site i. An external magnetic field will alter 
the entanglement between the carrier's spin and the LES, 
which leads to a magnetic-field modulation mechanism. 

A further issue needs to be addressed if the motion of 
the carrier's spin is treated as coherent, in which case 
the composite system evolves with time, and thus the 
CT rate, i.e., 77 = 77(f) from pi = pi(t). To reach a time- 
independent quantity for the effective rate, we treat the 
CT process as Markovian. Suppose the probability that 
the carrier is still on site i at time t is designated as 
P(t), the equation for P(t) is dP(t)/dt = -v r)(t)P(t). 
The average waiting time for the carrier on site i is 



(t w ) = 



exp 



r q{r)6.r I dt 



(4) 



If the contribution from the present mechanism is ne- 
glected, rj(t) = 1, Vq = l/(t w ). Similarly, an equivalent 
quantity v = l/(t w ) can be defined as the effective CT 
rate when 77(f) is incorporated. 

To present the theory quantitatively, we consider a ba- 
sic CT process, which is a charge carrier hopping from a 
site i to its neighboring vacant site j. This process could 
influence both the charge transport and the electrolumi- 
nescence. We take a hyperfine interaction model for the 
coupling between the carrier's spin and the LES, which 
is O 



H, 



(•5) 
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where we have set U = 1; a is the label for the differ- 
ent constituents of the LES, with J a the corresponding 
coupling strength and l a the nuclear spin operators; s 
are the carrier's spin operators; g is the Lande factor 
which is taken to be 2; pb is the Bohr magneton. The 
magnetic field B is chosen to be along the z-axis with- 
out loss of generality. The main result is shown in Fig. 
[U where the relative change of the hopping rate with 
the magnetic field is Au/v = \v{B) — j/(0)]/f(0). In the 
calculation, we have taken one spin-^ mode as the LES 
and set vq = 3.5 J. The final result is an average over 
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FIG. 1: The magnetic-field dependence of hopping rate Au/v 
with i/o = 3.5 J. The dashed line is the fit by a Lorentzian line- 
shape for the saturation component, with B1/2 the half-width 
at half-maximum (HWHM). The inset shows the ultrasmall- 
field component of Av/v, where B m is the magnetic field un- 
der which the hopping rate reaches its minimum value. 
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FIG. 2: The i/o-dependence of the characteristics of the 
magnetic- field dependent Avjv. The variation of the sat- 
uration value (Af/i/) S at and the minimum value (A^/V) m i n 
with vo are shown in (a); the variation of the HWHM B 1 / 2 
is shown in (b) and the variation of the magnetic field un- 
der which the hopping rate reaches its minimum value B m is 
shown in (c). The data of B m has been multiplied by 10 for 
clarity. 



all the initial states, with the carrier's spin taking any 
orientation and the LES in any of its eigenstates in the 
external magnetic field. Two components are present in 
the magnetic-field dependence of the effective hopping 
rate v. a saturation component at large magnetic fields 
which can be fitted well by a Lorentzian-like lineshape 
B 2 1 (B 2 + B 2 ^ 2 ), together with a negative ultrasmall-field 
component. This result is in good agreement with the 
experimental observations for MR [17j, as well as MEL 
[19{ . We emphasize that this kind of magnetic- field de- 
pendence is the general outcome of our theory, which 
does not depend on the choice of parameters. For exam- 
ple, we show in Fig. [2] the variation of the characteristics 
of this magnetic- field dependence with vq- Four quan- 
tities are calculated, including the saturation amplitude 
(Ai//u) max , the half- width at half-maximum (HWHM) 
B1/2, the minimum value of the ultrasmall-field com- 
ponent (A^/^) m ; n and the corresponding magnetic field 



B m . When i/ increases, the integrand of Eq. Q decays 
more quickly and there is less time to build up the en- 
tanglement, so both (Au/u) sat and [Av /v) m ^ n decreases. 
The ultrasmall-field component vanishes for Vq beyond 
about 5 J, because this component is from the oscillation 
behavior of the entanglement, which will not take effect 
for large vq. Furthermore, while B m remains nearly un- 
changed, B1/2 increases steadily with vq. We also point 
out that the two-component magnetic-field dependence 
of Avjv persists when more spin-i modes are included 
as the LES, only with quantitative modifications of the 
above quantities. The results for the case of two un- 
corrected spin— i modes as the LES are shown in the 
supplementary material [23|. 

The two-component magnetic-field dependence can be 
understood through some analytical results. We consider 
four basic initial states \a, 6), with a,b — ±i, in which a 
for the carrier's spin and b is for the LES. Other initial 
states can be seen as a proper superposition of the four, 
so the situation of the four could reflect qualitatively the 
behavior of Fig. [TJ from detailed averaging procedures. 
With the hyperfine-interaction Hamiltonian of Eq. ([3]), 
the total spin in the z direction, S z = I z +s z , is conserved, 
and we can classify the initial states by this quantity. 
Both the initial states of 3,5) (S z — 1) and I — |,— |) 
(S z = —1) are eigenstates without entanglement, giving 
i](t) = 1. However, the other two initial states — 5, 3) 
and 3,— 3) {S z — 0) become entangled states evolving 
with time, giving the same nontrivial result 

3 sin 2 cot + 2a 2 (I - coswt) 



where u) = Jyl + a 2 and a = g^sB/J. A time- 
independent quantity can be got by taking the time- 
average, which is 

, x 3 a 2 (a 2 - 2) 

5 (1 + a A ) z 

where we have scaled the quantity so that f\(B = 0) = 1. 
This direct time-average corresponds to the situation 
when Vq is sufficiently small, and the waiting time is much 
longer than the oscillation period in rj(t). Eq. ([7J is plot- 
ted in Fig. G2a), and it is clear that fj is less than 1 for 
a < y/2, which gives the negative ultrasmall-field com- 
ponent, and saturates for large B. This behavior can be 
traced back to the two-fold role of the increasing Zeeman 
splitting with the magnetic field on the entanglement. 
For example, for the initial state |^, — |), with magnetic 
field being zero, the initial state could evolve to the unen- 
tangled state \ — \, \) halfway in a complete period; with 
small magnetic fields, the degeneracy between the two 
untangled states is lifted, and the system remains entan- 
gled; with sufficiently large magnetic fields, the energy 
cost of flipping the carrier's spin allows no sizable prob- 
ability amplitude of | — 3, 5)- The entanglement remains 
small and the hopping ability becomes larger. 
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FIG. 3: (a) The magnetic-field dependence of the time- 
averaged fj(B) for each one of the four states (\a, b), a, b = ±4) 
as the initial state. The result of the Sz = states show 
both the negative ultrasmall-field component and the satu- 
ration component at large magnetic field; (b) The magnetic 
field dependence of 77(B) for the eigenstates of the hyperfine- 
interaction Hamiltonian, which only shows the saturation 
component at large magnetic field without the ultrasmall-field 
component. 



Till now, we have been considering the situation of the 
coherent spin motion on the initial site i. However, it 
should be noted that in some cases the coherency might 
be disturbed, such as when the spin-orbit interaction is 
enhanced by doping transition metal complex [24j. Al- 
though the interaction becomes more complicated, the 
results based on the hyperfine-interaction model of Eq. 
(O might still be indicative. In Fig. [3]Jb), the magnetic- 
field dependent fj(B) is shown for the eigenstates of the 
Hamiltonian of the composite system consisting of the 
carrier's spin and the LES, which are the states when 
the composite system reaches thermal equilibrium. Sim- 
ilar to the previous situation, two of the four eigenstates, 
which are ||, 5) and —5, —5) with S z = ±1, are unen- 
tangled states with rj(B) = 1. The remaining two eigen- 
states with S z — are entangled ones with 



fj(B) = 1 



3a 2 



(8) 



It can be seen that the ultrasmall-field component is not 
present, in contrast to the results for coherent spin mo- 
tion, while the saturation component survives the loss 
of spin coherence. This implies that the ultrasmall-field 
component observed in experiments could be a sign of the 
spin coherency. The persistence of the saturation compo- 
nent might also explain why the magnetic-field effect is 
still observed in the devices made by hydroxyquinolatcs 
consisting of heavy metal atoms [25] . 

We further calculate the isotope effect of the magnetic- 
field dependence of the hopping rate, by replacing the 
spin— i protons by the spin— 1 deuterons, and the result 
is shown Fig. |4] The coupling strength for the spin— 1 
case is chosen to be 0.5 J, where J is that of the spin— \ 
case, as the mass of deuterons are twice that of ordinary 
hydrogen atoms. The resulting HWHM of the spin— 1 
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FIG. 4: The isotope effect obtained by replacing the spin— | 
proton with the spin— 1 deuteron as the LES. The behavior 
at ultrasmall magnetic field is shown in the inset. Compared 
with the spin— j case, the ultrasmall-field component for the 
spin— 1 case becomes almost flat. 



case is larger than that of the spin— 1/2 case. Further- 
more, only a tiny ultrasmall-field component is present 
for the spin-1 case, which is shown more clearly in the 
inset. The reason behind the difference can be under- 
stood by similar analysis as above, which is presented 
in the supplementary material |23ij. Without any fur- 
ther assumption, the results obtained here is in excellent 
agreement with the isotope effect of MEL observed ex- 
perimentally [19j . 

Finally, we would like to emphasize that this work aims 
to present a general microscopic mechanism describing 
how the incoherent CT processes themselves are influ- 
enced by the carrier's spin states. With this mechanism, 
we can understand the peculiar properties of the MR 
and MEL observed for many different organic materi- 
als. Quantitative comparison with specific experimen- 
tal curves is beyond the scope of this paper and could 
be reached by considering specific materials as well as 
the working conditions. Besides the hyperfine-interaction 
model, this mechanism can be applied to other CT pro- 
cesses depending on which process is dominant in the 
specific organic devices. In some cases, there might exist 
more than one dominant CT process, when the magnetic- 
field effect should display multiple Lorentzian-shape com- 
ponents. This kind of behavior has actually been con- 
cluded by the careful analysis of experimental data [26j]. 

In summary, we have proven that the entanglement be- 
tween the carrier's spin and the local environment con- 
tributes a factor to the incoherent CT rate. This offers 
a general magnetic-field modulation mechanism, which 
leads to the MR in organic semiconductors. For the hop- 
ping of a single carrier to a vacant site, by considering a 
hyperfine-interaction model, both the saturation compo- 
nent and the negative ultrasmall-field component emerge 
naturally. The mechanism holds promise for the incorpo- 
ration of other influential CT processes leading towards a 
more comprehensive understanding of the magnetic-field 
effect in organic materials. 
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